Large-scale perturbations on the brane 
and the isotropy of the cosmological singularity 
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We present the complete set of propagation and constraint equations for the kinematic and non - 
local first order quantities which describe general linear inhomogeneous and anisotropic perturba- 
tions of a flat FRW braneworld with vanishing cosmological constant and decompose them in the 
standard way into their scalar, vector and tensor contributions. A detailed analysis of the perturba- 
tion dynamics is performed using dimensionless variables that are specially tailored for the different 
regimes of interest; namely, the low energy GR regime, the high energy regime and the dark energy 
regime. Tables are presented for the evolution of all the physical quantities, making it easy to do 
a detailed comparison of the past asymptotic behaviour of the perturbations of these models. We 
find results that exactly match those obtained in the analysis of the spatially inhomogeneous G2 
braneworld cosmologies presented recently; i.e., that isotropization towards the Tb model occurs for 
7 > 4/3. 
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I. INTRODUCTION 

A well known problem of cosmology is to explain the 
very high degree of isotropy observed in the Cosmic Mi- 
crowave Backgromid (CMB). In a theory such as gen- 
eral relativity, where isotropy is a special rather than 
generic feature of cosmological models, we need a dynam- 
ical mechanism able to produce isotropy. Inflation was 
proposed, among other reasons, as a way to isotropize 
the universe. Inflation is effective in this sense, but it 
needs homogeneous enough initial data in order for in- 
flation to begin 0|. Although one could perhaps adopt 
the view that one smooth enough patch in an otherwise 
non - smooth initial universe is all that is needed to ex- 
plain observations, this may not be satisfactory 0: the 
isotropy problem remains open in standard cosmology. 

Recently, a number of authors 0111111111111111 
ITlL lll | have addressed the issue of isotropization in the 
context of braneworld cosmology based on a generaliza- 
tion of the Randall and Sundrum model Il4j. Here 
the bulk is 5 - dimensional and contains only a cosmo- 
logical constant, assumed to be negative (see for a 
comprehensive review). 

In all cases considered, an interesting result was found: 
unlike general relativity, where in general the cosmolog- 
ical singularity is anisotropic, the past attractor for spa- 
tially homogeneous anisotropic models in the brane is a 
simple Robertson- Walker (RW) model |10|- Since 
this result was also found to hold for Bianchi IX models 
0,01 as well as for some inhomogeneous models, the au- 
thor suggested that the isotropic singularity could be a 
generic feature of brane cosmological models. 

In a recent paper 0, this conjecture was supported 
by studying the dynamics of a class of spatially inhomo- 



geneous G2 cosmological models in the braneworld sce- 
nario. A numerical analysis of the governing system of 
evolution equations led to the result that for 7 > 4/3 
isotropization towards a simple RW model Tb occurs as 
T —00 for all initial conditions. In the case of radia- 
tion (7 = 4/3), the models were still found to isotropize 
as r —>■ —00, albeit slowly. It can therefore be con- 
cluded that an initial isotropic singularity occurs in all of 
these G2 spatially inhomogeneous brane cosmologies for 
a range of parameter values which include the physically 
important cases of radiation and a scalar fleld source. 
The numerical results were confirmed by a qualitative 
dynamical analysis and a detailed calculation of the past 
asymptotic decay rates IT]. 

A similar result is also obtained in a related pertur- 
bative study where a careful analysis of generic linear 
inhomogeneous and anisotropic perturbations of the JT^ 
model [ill was conducted. Solutions were obtained for 
the large-scale evolution of scalar, vector and tensor per- 
turbations showing that the J-f, model is stable in the 
past (as r — > —00) with respect to generic inhomoge- 
neous and anisotropic perturbations provided the matter 
is described by a non - inflationary perfect fluid with 7 - 
law equation of state parameter satisfying 7 > 1. In 
particular, it was shown that the expansion normalised 
shear vanishes as r ^ —00, signalling isotropization. 

Brane cosmology thus has the very attractive feature 
of having isotropy built in, and although inflation in this 
context would still be the most likely way of producing 
the fluctuations seen in the CMB, there would be no need 
for special initial conditions for it to start. Also, the 
Penrose conjecture 19] on gravitational entropy and an 
initially vanishing measure of the Weyl tensor might be 
satisfied, c.f. 12(1. 
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The aim of this paper is to give a more compres- 
sive large - scale perturbative analysis of flat Friedmann - 
Robertson - Walker (FRW) brane models with vanishing 
cosmological constant by combining the high energy re- 
sults in _18J with an analysis of the other important stages 
in the braneworld evolution, namely the low energy GR 
and dark energy regimes. To make this precise we de- 
fine dimensionless variables that are specially tailored for 
each regime of interest. In this way we are able to clarify 
further the past asymptotic behaviour of these models 
and obtain results which match the analysis of the spa - 
tially inhomogeneous G2 cosmologies presented in |l7| : 
i.e., that isotropization towards the Th model occurs for 
7 > 4/3. 

The paper is organised as follows. In section II we 
will give a brief summary of the braneworld scenario and 
the induced field equations on the brane. In section III 
we introduce dimensionless expansion normalised vari- 
ables and derive the complete set of propagation and con- 
straint equations for the kinematic, inhomogeneity and 
non-local quantities. In section IV we split these equa- 
tions into scalar, vector and tensor parts, which we then 
analyse and discuss in sections V - VII for the low en- 
ergy, high - energy and dark radiation dominated regimes 
respectively. Finally, in section VIII we present our con- 
clusions. For the most part we follow the notation and 
convention of [isl IT^ . 



II. BRANE DYNAMICS 



equations above are not independent, but related to each 
other by 



A = 6- 



A= i[A(5) -Kk^A], 



(4) 



The tensor Sab represents non- linear matter corrections 
given by 



q J_7^^ T' — LT' T'^ -I L/^ 

'-Jab — 12 c-^ab A-'-ac-'- b i OA-J'^ 



cd 



(5) 

is the projection of the 5-D Weyl tensor C 



''ah 

on to the brane: E^^^ = C*-^^ 



-,C^D„A„B 



(5) 

(5) 

ABCD 

ABCD"- 9a9b^ where ua is 
the normal to the hypersurface x = {uAn^ = ^1)- 
Although the whole dynamics are 5-D and given by 
, from the 4 - D point of view i?^^-* is a non - local source 
term that carries bulk effects onto the brane. 

The energy - momentum tensor Tab is assumed to be 
conserved on the brane: 



V'Tab = , 



(6) 



and on using the 4-D contracted Bianchi identities 
V'Gab = an additional constraint is obtained: 



(7) 



which shows how the non - local bulk effects are sourced 
by the evolution and spatial inhomogeneity of the brane 
matter content. 



A. Geometric Formulation 

The implementation of the braneworld scenario consid- 
ered in [ij] assumes that the whole spacetime is 5 - D and 
governed by the 5-D field equations {A,B = 0, ...,4): 

g'^ab = -A(5).9!fi + «(5)'5(x)[-AgAB + Tab] • (1) 
These represent a 4 - D brane at x = embedded in a 

(5) 

vacuum bulk with metric g^^ and cosmological constant 

A(5); 

is the 5-D gravitational constant, A is the 
brane tension, gAB and Tab are respectively the met- 
ric and the energy- momentum on the brane. The 4-D 
field equations induced on the brane are derived geomet- 
rically from ^ assuming a Z2 symmetry with the brane 
at the fixed point, leading to modified Einstein equations 
with new terms representing bulk effects: 

Gab = -Agab + n^Tlf , (2) 

where 

T'^' = Tab + f - ^E'i^ . (3) 

As usual = Sn/Mp, and {a,b = 0, ...,3). The vari- 
ous physical constants and parameters appearing in the 



B. Cosmological Dynamics on the Brane 

In the following we describe the matter on the brane by 
a perfect fluid with barotropic equation of state p — {'y — 
l)p. As usual, we require 7 > to satisfy the dominant 
energy condition and 7 < 2 to preserve causality and 
therefore we restrict our analysis to values of < 7 < 2. 
The case 7 = can be treated similarly to our analysis 
below, but using different variables. We do not study this 
special case here, but refer the reader to [23| for details 
on how to treat this case. 

If is the matter 4 -velocity and hab = gab + ""a^b 
projects into the comoving rest space of a fundamental 
observer, the brane energy momentum tensor is given by 

Tab^ pUaUb+phab ■ (8) 

One can also decompose E^f^ in such a way that it is 
equivalent to a trace-less energy momentum tensor with 
energy density p* , energy flux q* and anisotropic pressure 
TTab (see for details): 

- I^E'S' = P*{UaUb + \hab) + q*aUb + g>a + Kb ■ (9) 

Since there is no evolution equation for the non-local 
anisotropic pressure 7r*j, we restrict our analysis to the 
case 7r*j, = 0. This is dynamically justified in early time 
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regimes Note in particular, that using (O 

the total energy density is given by 



(10) 



In this way we can see that there are three essentially 

2 

different energy regimes: when p ^ \^P* recover 

2 

general relativity (GR). When ^ '3> p,p* we obtain the 



high energy limit, and when p* 3> p, ^ we obtain the 



dark radiation dominated regime. 

Equation 10 gives the usual energy and momentum 
conservation equations: 



p + 3jHp = , 

(7- l)D„p + 7pAa = 0, 



(11) 
(12) 



where a dot denotes u''Vb, H = SD^Uq is the Hub- 
ble parameter of the background, Ai, = ut is the 4- 
acceleration, and Da denotes the spatially projected co- 
variant derivative. 

Using equations @ and (jZj) we can obtain conserva- 
tion equations for the non - local quantities p* and q* . 
Restricting to linear perturbations of Robertson - Walker 
models illMilEi, we obtain Issll: 



P 



4iJp*-fD^: = 0, 



(13) 



ql + AHql + \^aP* + ip*Aa = -f 7DaP . (14) 

Finally we note that the generalised Friedmann equa- 
tion on the brane for a flat, homogeneous isotropic back- 
ground with vanishing 4-D cosmological constant A is 



3 6A 3 



(15) 



III. COSMOLOGICAL PERTURBATIONS 

In the following sections we will present a complete 
description of general large - scale inhomogeneous pertur- 
bations for several different flat {K = 0) homogeneous 
isotropic background models with vanishing cosmologi- 
cal constant on the brane (A = 0). 

Following the dynamical systems approach developed 
in I3ljj32l 1331 and extended to the braneworld scenario 
in IJ, ly] we define the dimensionless density parameter 
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3IP 



as in general relativity, and 



p^ hl^ 



3iJ2 



(16) 



(17) 



corresponding to the non - GR contributions to the Fried- 
mann equation. In this way we can classify the various 
background solutions by their coordinates (p,pVLx,Q,p*) 
in the phase space of Friedmann- Robertson- Walker 
(FRW) models. 



The point (1,0,0) corresponds to the flat GR Fried- 
mann model with A'^ = p* = and a{t) ~ ^2/37^ The 
point (0, 1,0) corresponds to the high energy model Tb 
with p* — Q and p <^ p^/A; the scale factor is given by 
a{t) 'y t^/'^'*', which can be found by a limiting process 
Finally the point (0, 0, 1) corresponds to a model 
(i?) with scale factor evolution a{t) ~ i^/^, satisfying 
p = A-i = 0. 

In what follows, we develop the perturbation equations 
for a general background model {Vlp f^A, Op. ) and decom- 
pose the propagation and constraint equations into their 
respective scalar, vector and tensor contributions in the 
usual way. We then evaluate these equations for the three 
backgrounds described above and interpret the results. 



A. Dimensionless Variables 

The projected 4-D field equations 10) can be covari- 
antly S£lit using the Ricci identities and the Bianchi iden- 
tities |l5| . In the previous section we have already given 
the conservation equations for energy and momentum 
(|ll|l , (|12|l and for the non - local energy density p* and 
flux (7* ifT^ . The remaining equations correspond to 
propagation and constraint equations for the kinematic 
quantities i.e., the acceleration A^, the vorticity cjf, and 
the shear Gab, together with the electric and magnetic 
parts of the Weyl tensor Eab, Hab corresponding to the 
non-local gravitational field on the brane. 

Instead of using the standard quantities we define di- 
mensionless expansion normalised variables by 



Ti/ _ „ _ CTafc p _ ^ab nj _ ^ab 

= , 2^ab = — , Oab = , rLab = ■ (ioj 



H 



We emphasise that Eab and the expansion normalised 

quantity Sab niust not be confused with the 5 - D Weyl 

(5) 

tensor E^^^ . It turns out that using the dimensionless 
vorticity variable 



W: = aHWa 



(19) 



simplifies the calculations below; however, it should be 
noted that Wa and not W* is the physically relevant 
quantity. We also use the dimensionless logarithmic time 
derivative r, defined by 



d ln{a) 



Hdt 



(20) 



Density perturbations are physically characterised by 
the comoving fractional density gradient defined by |22l | : 



-DaP. 

P 



(21) 



In addition, it is convenient to define the following dimen- 
sionless gradients describing inhomogeneity in the expan- 
sion rate H and in the non - local energy density p* and 
flux ql 

z: ^ ^DaH , u: EE ^Bap* , q: ^ . (22) 
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Note that although Aa is not defined for the exact dark 
radiation background (i?) where p = 0, it is is well- 
defined in a neighbourhood of (i?) and since (i?) is a 
saddle point in the phase space of the background ho- 
mogeneous models Jj] , this solution can never be exactly 
attained. It is therefore sufficient to study arbitrarily 
small but non - zero inhomogeneous perturbations of this 
model. In other words we may evaluate the perturbation 
equations arbitrarily close to the background (i?), but 
not on the exact background itself. 

The above discussion suggests that it makes more sense 
to define specially tailored inhomogeneity variables by 
normalising them with respect to the dominant energy 
density term in the Friedmann equation. 

For the low - energy limit we use the usual dimension- 
less perturbation quantities : 

A^) ^ ^D,p, C/r) ^ ^D^p*, Q^) ^ , (23) 
P P P 

however since is the dominant term in the dark radi- 
ation dominated regime, the appropriate inhomogeneity 
variables are 

A(de) ^ o U^^^^ = — D o* O^^^^ = —a* 



Finally in the high - energy limit we define 
Aa 



(24) 



p 



-^BaP 



P 



(HE) 



P^ 



-7^1^ ■ 
(25) 

since the leading energy density term is proportional to 
p^ l\ which becomes our normalisation factor in this 
case. 

When decomposing the equations into harmonics we 
will have to deal with the curls of some of the variables. 
One approach would be to eliminate the curls by deriving 
second and higher order equations, or alternatively intro- 
duce new spatial harmonics corresponding to the curls of 
the original harmonics. Instead, we find it more conve- 
nient to define new variables corresponding to the curls 
of the original quantities and derive propagation and con- 
straint equations for them. In this way we obtain a com- 
plete closed set of linear differential equations which can 
be easily solved. Note that all of the additional propaga- 
tion and constraint equations have to be satisfied, since 
these equations are necessary to close the system. 

We denote the curl of a quantity with an overbar and 
the key variables of this type are: 



and 



Wl = -|curW; , = -^curlS.fa , 

^ab = -yrcml Sab ,'Hab = -^cmlHab 
Id Id 



Q:^-curig:. 



(26) 



(27) 



B. Dimensionless Linearized Propagation and 
Constraint Equations 

The complete set of propagation and constraint equa- 
tions for the kinematic and non-local quantities on the 
brane were developed in [l5| . Here we extend this work 
by presenting the complete set of evolution and con- 
straint equations for the dimensionless variables defined 
by frmm and (PI). 

We begin with the generalised Raychaudhuri equation 

H' , 

(1 + 9) 



7-1 



D"A„ 



(28) 



where 



q = |7l7p + 37f^A + 2f7p. - 1 



is the usual deceleration parameter. 

The remaining propagation equations are given by 

w:' - (37 - 4)1^: , 

= (2g - l)£ab -iq + l)Sa6 

~ 2aH'^ D<aQ*h> + T^ab , 
n'ab = (2<7 - l)nab - £ab , (29) 

which are subject to following dimensionless constraints: 



oD^w; = , 



2 7* 
3^a 



a-D'Hab = 2(q + 1)1^: - ig: , 

where the angle bracket is defined by 



\D'W:hab 



(30) 



(31) 



Equations for the inhomogeneity variables A^, Z*, Q* 
and Q*. are given by 



A,' 


= (37 


- 3)A, - 7Z: 










i)z:-|[i7p-f 


-(67+2)r!A- 






-u:- 


6(7 - l)Wa - 






Qa 


= iq- 


2)q: - ^,u: 4 


- (4^f^p. - 





mp,zi 



u:' = {2q - 2)u: 

+ i^rip- A, - ^D,{B'Ql) 



(32) 



Finally, using the definitions H2t)|l and (|27|l . we obtain 
equations for the curls of the original quantities: 

W*'a = (9 + 37 - 4)14^: , 

s;., = {2q - i)s,b - Bab , 

E'^^ = (Sq - l)£ab - (g + l)Sa6 + ^r>^aT^'Hb>c 
^Ib = (39 - ^)'Hab - 2^D<qD^'^6>c + -jJV^^Sab , 

g*l = (2q - 2)Q: + 4[(67 - 8)17p. - 9^^nx]W: (33) 
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which are subject to the following constraints (obtained 
by taking the curls of (|3U|l '): 



Z: = (q-l)Z,-[|f]p + (97 + 3)f7 



(39) 



aWw;; = , 

aD^S,, = 2{q + 1)W: - \Ql - ^V>^W: , 



aD'£ab^2iq + l)w:-^Q:, 

aD'Hab = {q + l)W: + 4^(D2q: - D,(D''Q^)) , 

aD''Q: = 0. (34) 

IV. HARMONIC DECOMPOSITION 

In order to solve these equation we employ the stan- 
dard approach of expanding the variables in these equa- 
tions in terms of scalar (S), vector (V) and tensor (T) 
harmonics Q js^ . These harmonics are eigenfunctions of 
the covariantly defined Laplace - Beltrami operator : 



(35) 



where k is the wave number corresponding to a comov- 
ing scale A = 2'Ka/k. This yields a covariant and gauge 
invariant splitting into three sets of evolution and con- 
straint equations for scalar, vector and tensor modes. 

Thus a scalar X, vector X° (orthogonal to u") and 
tensor X""^ (orthogonal to u") can be expanded as follows 



X 

Xa 

Xab = fc^'^^Qfb 



k-'X^'Q^. + X^Ql. (36) 



In what follows we drop the subscripts 5, V, T and also 
restrict our analysis to the long wavelength limit defined 



by 



« 1. 



Scalar Perturbations 



In the long wavelength limit the scalar evolution equa- 
tions for the kinematics (which follow after expanding 
(|29|l in terms of scalar harmonics) are given by 



E' = (g - 1)S - £ , 

= (2g - 1)5 - (g + 1)S 



(37) 



and these are subject to the following constraints (which 
follow from 

3(f^p + 2f7A)A = 2£: + 3Q, - . (38) 

In addition, we have scalar evolution equations for the in- 
homogeneity variables (which follow from (|32|) after har- 
monic analysis): 

A' = (37 - 3)A - , 



6^r!p. ] A - [/, 



= (g - 2)Q, - \U, + {A^i^^p. - 6717a) a , 



U', = (2q - 2)C/, - AVLp, + U^Vtp, A . (40) 



B. Vector Perturbations 

Expanding equations (|29l) and (|33|l in terms of vector 
harmonics, we obtain the following evolution equations 
for the kinematic and non - local quantities together with 
their curls: 

Wj = (37 - 4)1^, , 
E' = (g - 1)S - 5 - 2(7 - 1)T^, , 
£' = {2q-l)£ - {q+l)J: + H , 
W = {2q-l)n-£ , 

= (37 - 4 + q)W, , 
E' = {2q-l)f:~-£ , 
£' = {3q-l)£-{q+l)E , 

H' ^ i3q-l)H , 

= {2q - 2)Q, + 4[(67 - 8)17p- - 97^^!^]^^* • 

These equations are subject to the following constraints, 
which are obtained from H30|l and (|34|l : 

= 3S] -f 6Q* - 6W, , 
6{np + 2nx)A = 3£ + 6Q, - 2U, , 
S = H , 

n^i{q+l)W,-Q, 
. H^2{q+l)W, = , 
E = 4(g + 1)W, - Q, , 

£ = A{q+\)W.-Q, . (41) 

Equations for the vector parts of the inhomogeneity vari- 
ables follow from (|32|l and are given by 

A' = (37 - 3)A - 7Z, , 

Zj = {q- 1)Z, - [^Qp + 3(37 + l)f^A - 6^f]p. ] A 
-C/, - 6(7 - 1)W, , 

QJ = {q- 2)0, - if/, + (42^r!p. - 67S1a)A , 

= {2q - 2)U^ - mp, Z^ + 12^f7p. A . (42) 



C. Tensor Perturbations 

Finally the long wavelength behaviour of tensor per- 
turbations are obtained by expanding H29|) and (|33f) in 
terms of tensor harmonics: 

E' = (g-l)E-£, 

5' = (2g-l)£:-(g+l)S, 
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t' = {2q-l)t,-£ , 

£' = {3q-l)£ - {q + l)t , (43) 
subject to the following constraints 

7i = , E = H . (44) 



V. LOW ENERGY LIMIT: THE GR 
BACKGROUND 

We begin with perturbations in the low -energy limit, 
defined by p P^/A and p :s> p* or ftp V,x and 
rip ^ Hp* . We therefore evaluate the perturbation equa- 
tions H38I43|I in the limit {np,nx,Hp') ^ (1,0,0). Us- 
ing the energy conservation equation the Friedmann 
equation H15|l can be solved to give the background scale 
factor a and the Hubble parameter H: 



a{t) = (t/to)'^^^ , H = Hoa~^^'^ 



(45) 



where we fix an arbitrary initial condition by choosing 
flo = a(to) = 1- The deceleration parameter is given by 



and, as usual, 



rj 2 

where = ^po- 



9= 47-1, 



p = poa ^'^ , 



Scalar Perturbations 



(46) 
(47) 



In this case the propagation equations for the kine- 
matic quantities are 



(§7 



2)S-£. 



£' = (37 - 3)£ - I7E , (48) 
while the constraints are given by 

w = w = H = n^f:^£^Q^o, 

2E + 3Q* , 



3A = 25 + 3Q* - [/* 



(49) 



and the equations for the inhomogeneity variables be- 
come 



A' = (37 - 3)A - 7Z, , 

= (§7 - 2)Z, - |A - C/, , 
g,' = (|7-3)Q, - iC/*, 
Uj = (37-4)C/, . 

The above equations can be easily solved to give 



E = ^oa^'^-^ + Eia2 



7-3 



j7Eoa 



37-2 



+ Sia2 



7-3 



and 

A = —7^00' 



37-2 



(IE, 



Eoa3^-2 



Ql)a2 
3 



+ (Ei + fQ5)a2 



7-3 
7-3 



7 rr*„37-4 



37-2 



37-2*^0 



2 

3 37-2 



3— r/*o^'''~'' 



(52) 



for 7^1, and 



A = -|Eo + (|Ei+i 



= Eo + (El + |gS)a-2 - iC/o*lna a"^ ^ 



i C/q In a a 



(53) 



if7 = f. 

Here Eq, Ei, Qg, Uq are arbitrary constants of inte- 
gration, corresponding to the four independent modes. 

Finally the solutions for Q*, [/* can be converted into 
the the physical quantities Q^^^\ [/(le) (-^^iijch corre- 
spond to the scalar modes of Qa, Ua defined in (|23() ) 
using the background solutions for H and p: 



Q 



(LE) 



LE) 



irr*„37-4 
3^0" 



for 7 7^ I and 



g(LE) _ ^^^g*Q 
{/(LE) ^ 1 [7*^-2 



9ffo 37-2*^0" 



(54) 



(55) 



if 7 



B. Vector Perturbations 



For vector perturbations the complete set of propaga- 
tion equations for the kinematic and non - local quantities 
are given by the ten - dimensional system 

Wj = (37 - 4)W, , 
E' = (|7-2)E-£, 
£' = (37 - 3)£ ~ I7E , 



n' = (37 -3)n~£ , 

E' = (37 - 3)E - f , 
£' = (§7 - 4)5 - f 7S , 
Q'. = (37 - 4)Q* . 

(50) subject to the following constraints 

4z, = 3E + eg* , 
6A = 3f + eg, - 2c/, , 

(51) n = w, = o, 



(56) 



(57) 
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TABLE I: Large scale contributions of the different modes to the geometric and kinematic quantities for the low energy 

background. We assume 0<7<2, 77^2, and wc omit non-zero constant coefficients. The first line in this tabic should be 

read as E = aSoa^^~^ + /3Sia2^~^, where a, (5 are some non-zero constants, for all scalar, vector and tensor modes. General 
Relativity is recovered when Qq = = 0. 



harmonic 


scalar 


vector 


tensor 


mode 


a'-'-^ ah-^ a'-f-^ ah-^ 


^3^-2 ^|t,-3 ^37-4 ^§7-5 


^37-2 af7-3 ^37-4 ^|t-3 


S 
£ 

n 
w 

A 

z* 

g(LE) 
jy(LB) 


Eo El - - 
Eo El - - 

Eo 2S1+3QS US 
Eo 2S1+3QS US 

Ql US 

US 


Eo El - - 
Eo El - - 

ws 

Eo E1+2QS US 
Eo S1+2QS US 

Ql US 

US 


Eo El - 
Eo El - 

Tio "Hi 



while the propagation equations for the inhomogeneity 
variables are 

A' = (37 - 3)A - 7Z, , 
ZJ = (§7 - 2)Z, - |A - [/, , 
Qj = (§7-3)0, - \U, , 
Uj = (37 - 4)C/, . (58) 
Solutions can again be easily obtained and are given by 



37-2 



7-3 



= E = f = Hoa 



37-4 



(67VK0* - Ho)a 



37-4 



(59) 



and 
A 



^ — 3Ho^O"- 9Ho 37-2^0" : 

?7(LE) = l[/*a37-4 



(60) 



for 7 7^ |, and 



^Uq In a a 



Z, = fSo + (fSi + |0S)«"' - Ina a'^ 



(61) 

if 7 = |. Again, we have converted the solutions for 
Q^, into the physical quantities Q^^^\ U^^^K 

There are six independent modes corresponding to the 
constants of integration Sq, Si, Wq, Hq, Qq and Uq. 



C. Tensor Perturbations 

For tensor perturbations the propagation equations in 
the long wavelength limit are 





= ih 


-2)S- 


-s, 


£' 


= (37- 


-3)£- 


Its, 


n' 


= (37 


-3)n- 


-s, 




= (37 


- 3)S - 


£, 


£' 


= (§7 


-4)£- 


i7S, 



subject to the following constraints: 
H = 0, S = W. 

The solutions are 



(62) 



(63) 



S = Eoa3T-2-HSia2''-3, 
£ = -|7Eoa3T-' + Elai^-^ 
H = f, = Hoa^^-^ + Wia2T-3 , 



(64) 



where Sq, Si, Hq, Wi are four independent constants 
of integration. 



VI. THE DARK ENERGY ERA 



The dark energy dominated regime is characterised by 

p* > p and p* » p^/A or flp, > Hp and Hp* fix, so 
we now evaluate the perturbation equations in the limit 

(fip,fiA,f^pO^ (0,0,1). 

The background solution (i?), has the same metric as 
a flat radiation FRW model, with p = and p* = PQa~^. 
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The scale factor a, Hubble parameter H and deceleration 
parameter q are given by 

a{t) ^ [t/tof^^ , H = Hoa-^, q = l. (65) 

As explained in section III the perturbation equations 
are only defined for small but non - zero energy density 
p, or equivalently arbitrarily close but not on the exact 
background (i?). We therefore use p = poa~^^ for small 
but non - zero values of po ■ 



A. Scalar Perturbations 

In the case of scalar perturbations the propagation 
equations reduce to: 



f = f - 2S , 
subject to the following constraints: 



(66) 
(67) 



w^w^n^n^Ts = £ = Q = Q, 

2£ + 3g, - [/, = 3(17p + 2f^A)A , (68) 

while the equations for the inhomogeneity variables are 

A' = (37 - 3)A - , 
Zj = 6^A - , 

Q,' = -g*-if/*+4^A, 



The following solutions can then be obtained: 

E = Ega^ + Eifl""'" , 
£ = -2Eoa2 + Eifl-i , 



(69) 



(70) 



and 



A = Ao + Aia37-5 + ^Sga^ , 
= 22_3 Ao + f Aia3^-5 + 3f^Eo 



7 

7 



37 



Aia^^-^ + lf^Eoa^-fEia-i, 

(71) 



[/. = ^Ao + |Aia3^-5 + gi^Eo 



for 7 |. We do not give the solutions for 7 = |, since 
all values of 7 > 2 are outside the region of interest. 

There are four constants of integration Eo, Ei, Aq, Ai 
corresponding to the four independent modes. 

The scalar contributions to the physical quantities de- 
fined in H24|l can then easily be obtained: 



= I^A 



2 



A ^37-4 I 4 37-4 ^ : 



(72) 



In particular, the density perturbation A^^^^ can be 
written as fJp. A*^^^^ = f^pA, and using the fact that 
Q,p,Vlx > 0, we find that ilpA — flxA — > as p — > 0. 
Hence the density perturbations arc given by 

A(DE) = ripAo + r!pAia37-5 + _|7_f^^Soa2 (73) 

and are suppressed as one approaches the dark energy 
solution. 



B. Vector Perturbations 

In the case of vector perturbations the propagation 
equations are 

= (37 - 4)iy, , (74) 

E' = -£, (75) 

£:' = £-2E, (76) 

H' = 1:' = £' = 0, (77) 

= 8(37 - i)W, , (78) 

which are subject to the following constraints 

3E + 6Q, = AZ, , 

35 + 60* = 2U^ + 6{Qp + 2nx)A , 

E = W, 

■w = 8iy* - 0* = , 
W = If, = , 

f: = £ ^8W,~Q, . (79) 
The inhomogeneity variables evolve according to 

A' = (37 - 3)A - 7Z, , 

= 6^ A - U. . 

1 

Q*' = -0*-iC/,+4^A, 

Uj = -4Z, + 12^A . (80) 

The solutions to this system are 

E = Eqo^ + Eia~^ , 
£ = -2Eoa2 + Eifl-i , 
7i — ^ — £ — Tig , 
W = 2Wo*a3'^-3 , 

Q, = 8W^a^^-* + Hn , (81) 

and 

A = Ao + Aia^^-s + fgl^Eoa^ , 

= i^Ao + f Aia^^-^ + ||25f Eoa^ 

= ^Ao + ^Aia37-5 + |25|E„a2 _ 1 ^^^-1 ^ 

(7. = ^Ao + f Aia3^-^ + 3^Eoa2 , (82) 
where again 77^^. 
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This time there are six constants of integration: where again we fix an arbitrary initial condition by choos- 

Sqi T^Qi Wq, Qq, Ai. ing ap — a{to) = 1. The energy density behaves in the 

The vector modes of the physical quantities defined in usual way: 
(I24|l can then be found: 

P = Poa'^'^ ■ (88) 

A (DE) _ K Po A „4-37 I K po A „-l I K Po 37 V „6-37 

^(DE) 27-2 A ,4a 37-4 , 4 37-4 V. 3 ^rom the Friedmann equation ^ we find that = 



C/(°E) = C/* . (83) 



Scalar Perturbations 



(84) 
(85) 



C. Tensor Perturbations 

The tensor parts of the propagation equations in the 
long wavelength limit are: 

f ' = £ - 21] , 

E' = S - £ , 
f ' = 25 - 2S , 

subject to the following constraints: 

71 = 0, t = n. 

The solutions are 

S = Soo^ + Sia""'" , 

£ = -2Soa2 + Sia"^ , 

Ti = Yi = Tio + Tiia^ , 

£ = no- 2Hia^ , (86) 

with a constant of integration for each of the independent 
modes: Eq, Ei, Hq, Hi. 



VII. HIGH ENERGY LIMIT: THE n 
BACKGROUND 



The high energy limit is characterised by ^ 



The scalar propagation equations for this case reduce 



and x > 



or Qx '3> flp and Qx ^ ^p' 



> P 
so this 

time we evaluate the perturbation equations in the limit 

{np,nx,np.) ^ (0,1,0). 

This model corresponds to a stationary (equilibrium) 
point ^ in the phase space of homogeneous Bianchi mod- 
els 0,0, as well as in the phase space of the special class 
of inhomogeneous G2 cosmological models. In both cases 
J-f, is found to be the source, or past attractor, for the 
generic dynamics for 7 > 1 (7 = 1 is also included in the 
homogeneous case) , consistent with Q |^ |^ . The sta- 
bility of this result is now examined through an analysis 
of the perturbation equations for this case. 

The background scale factor a, Hubble function H and 
deceleration parameter q of these models are given by 

a{t) = {t/tof'' , H = Hoa-^\ g = 87 - 1 , (87) 



to 



S' = (37 - 2)E - 5 , 
£' = (67 - 3)£ - 37S , 



(89) 



subject to the constraints 



W = W = H = H = T. = £ = Q = 0, 

6A = 2£ + 3Q, - . (90) 

The scalar evolution equations for the inhomogeneity 
variables are 

A' = (37 - 3)A - 7Z, , 

Zj = (37 - 2)Z, - 3(37 + 1)A - , 

QJ = (37 - 3)Q» - \U, - 67A , 

Uj = (67 - 4)C/, . (91) 

Solutions can again be easily obtained by solving the 
above system of linear equations. They are 



£ = -37Soa'^'^-2 _^ £^a37-3 



(92) 



and 



A - in*„-3 7(37+1) ^ 67-2 -y 7-r*„67-4 

^—2*^0" 67+1 ~ 2(67-1) "-^O" ; 

^* — 2^0" + 67+1 +2(67-1)^0" ' 

I J37-]_rr* 67-4 
^ 3(67-1)^0" : 

[/, = U^a^^-'^ (93) 
for 7 ^ ^, and 



A = (ig* - iC/o*)a-3 _ lEga-i - ^U* In a a'^^A) 
= |Oo«"^ + fSoa"^ - if/o In a a"^ , (95) 

Q, = QSa"^ + T2^oa-^ ~ fEia^t - i?7* In a a^i^) 
C/* = U^a-^ (97) 
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for 7 = g- The scalar parts of the physical quantities 
^(HE)^ q(he)^ jj(he) defined by ^ are given by 

A (HE) _ K^po n*„37-3 7(37+1) K^Poy „9f~2 

^ - rm^^oo- 67+1 617^^0°^ 

_ 7 K^po rr*r,97-4 
67-lT2i?J^0" ' 

I 1 37-1 rr* 97-5 
^ 18H0 67-1^0" 1 

[/(HE) ^ l[/*a67-4 (gg) 

for 77^^, and 

^(HE 



.(HE) 



1 



1 V 



(99) 



for 7 = i. 

So, Si, Qq, C/q are arbitrary constants of integration 
corresponding to the four independent modes. 



B. Vector Perturbations 

For vector perturbations the propagation equations are 

= (37 - 'i)W, , 
S' = (37 - 2)S - £ , 
£' = (67 - 3)£ - 37E , 
n' = (67 - 3)H - £ , 
S' = (67 - 3)S - f , 
f ' = (97 - 4)f - 37S , 

= (67 - 4)0, - , (100) 

subject to the following constraints 

= 3s + eg, , 

12A = 3£ + eg, - , 

n = t^£^ i2jw, - g, , 

n = w^ = o. (101) 

The propagation equations characterising the inhomo- 
geneities are 

A' = (37 - 3)A - 7Z, , 

= (37 - 2)Z, - 3(37 + 1)A - [/, , 

g,' = (37 - 3)g, - if/, - e7A , 

= (e7 - 4)C/, . (102) 



Solutions are easily obtained and are given by 

H = S = £ = 7^00^7-4 ^ 

= 37^0*0^'^"^ , 
g, = 127^0*03-^-" - Hoa^^-'' , (103) 



and 



a(HE) _ K.'^po ^, 37-3 7(37+1) K^Po V „97-2 

^ ~ mTJ^^oO 67+1 9i?^^"" 

_ 7 K^po rr*„97-4 
67~lT2i?J^0" ' 

n*„-3 I 3 (37+1)^^ ^ ,,67-2 , 37-1 rr* 67- 
'^ofi +4 67+1 ^0°^ + 2(67-1) "^o" 



/o(HE) _ ^r)*«3T-4 _L 1 37^ y,,r;9^-3 



4ffo 67+1 

1 : 

18^0 67-1 "^0' 



[/(HE) ^ i[/*a67-4 



(104) 



for 7 7^ i, and 



= |gSa-3 + §Soa-i - iWlna 0"=^ 



32' 



g(™) 

[/(HE) 



6Hq '3o" 



1 1 V 
677^T6^oa 



12_ff| 



■Eia 



-3 



if/*r)~3 



(105) 

(loe) 



for 7 = ^. There are six constants of integration: 
^0, "^07 Wq, Qq, Uq. 



C. Tensor Perturbations 

The tensor parts of the propagation equations in the 
long wavelength limit are: 





= (37- 


2)S- 




£' 


= (67- 


3)£:- 


37S, 


n' 


= (67- 


-i)n- 




t' 


= (67- 


3)S - 




£' 


= (97- 


4)£- 


37S, 



(107) 

subject to the constraints 

7i = 0, t^n. (108) 
The solutions are given by 

S Eoa^''"^ + Sia^'^-^ , 
£ = -37l]oa^^"^ + Sia3-y-3 , 
7^ = S = Woa^^"^ + ^la^^"^ , 
£ = Tioa®^""^ - iirLia^'<~^ . (109) 
There are four constants of integration: Eq: TYq, Hi. 
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TABLE II: Large-scale behaviour of the non-zero physicaUy relevant geometric and kinematic quantities for the different 
backgrounds in the limit a ^ 0. The values S, V, T in brackets denote scalar, vector and tensor contributions. One can easily 



see how the appropriately normalised perturbation quantities defined in l|23^-l|24|l 


converge 


for wider ranges 


of 7 as a 


0. 


Quantity 


mode 


0<7< i 


<'y<i 


< 7 < 




< 7 < i 


<7< 1 


<7<1 7<ia 


<7< 1 


< 7 < 1 


< 7 < 2 


Low energy limit 


E, £ 


^oiS,V,T) 


oo 


oo 


oo 


oo 


oo 






















Ei{S,V,T) 


CO 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


n 


no{V,T) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 










ni(T) 


oo 


oo 


oo 


oo 


oo 




















w 


w:;{v) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 











A = A(LE) 


Eo(S', V) 


oo 


oo 


oo 


oo 


oo 






















^i,Ql{S,V) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 




U*o{S,V) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 








Q(LE) 


Ql{s,v) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 












oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 











f/{LE) 




oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 








Dark energy limit 


E, £ 


MS,V,T) 





































Ei{S,V,T) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


n 


■Ho{V,T) 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 




Hi(T) 



































w 


w,;iv) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 














^(DE) a 


PoEo(5',V) 





































PoAo{S,V) 





























oo 


oo 




PoAi(S,y) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


g(DE) 


MS,V) 





































^i{S,V) 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 




Ao{S,V) 





































Ai(S,T/) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 








jj(DE) 


Eo(5', V) 





































Ao{S,V) 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 


const 




Ai{S,V) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 





High energy limit 


E, £ 


MS,V,T) 


oo 


oo 































Ei(S',\/,T) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 














n 


no{V,T) 


oo 


oo 


oo 


oo 


oo 






















ni{T) 


oo 


oo 





























w 




oo 


oo 


oo 


oo 


oo 


oo 

















/^(HE) 


Eo(5',V) 


oo 


































Qo{S,V) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 
















U5{s,v) 


oo 


oo 


oo 


























g(HE) 


Eo(5', V) 


oo 


oo 































El (5,1/) 


oo 


oo 


oo 


oo 


oo 






















QUS,v) 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 


oo 










us{s,v) 


oo 


oo 


oo 


oo 























[/(HE) 


U5{s,v) 


oo 


oo 


oo 


oo 


oo 





















"This quantity is suppressed by a factor of po, hence not significant 
when approaching the vacuum model (R). 
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VIII. RESULTS AND DISCUSSION 

In the previous three sections we have developed and 
solved the perturbation equations for the low energy, dark 
energy and high energy backgrounds, respectively. 

The main results of our analysis is summarised in Ta- 
ble II, in which we present the early time asymptotics 
a ^ of the physically relevant quantities for the dif- 
ferent energy regimes. These physically relevant quanti- 
ties are the harmonically decomposed components of the 
expansion normalised vorticity, shear, and the electric 
and magnetic parts of the Weyl tensor H18I19(I . as well 
as the appropriate gradients of the energy density p, the 
non - local energy density p* and the non - local flux q* 
defined in (|23(l - (|25|l . The remaining quantities appearing 
in the previous sections are required to close the system 
of equations, but are otherwise of no particular physical 
importance. 

We can see from Table II that in the low - energy regime 
the results from general relativity are recovered. In par- 
ticular, we find the same decaying mode Si in both the 
shear E and and density gradient A, implying that in 
general relativity the fiat RW models are unstable with 
respect to generic linear homogeneous and anisotropic 
perturbations into the past which was the problem out- 
lined in the introduction and partly the motivation for 
infiation. 

In the dark energy limit, we find that for any value 
of 7 there is a quantity that diverges as a ^ 0. The 
dark energy background is, however, an unstable equi- 
librium point in the state space of flat FRW models [4| , 



and can therefore only be attained for very special initial 
conditions. 

The main result of this analysis relates to the evolution 
of the perturbation quantities in the high energy back- 
ground JFf,. We find that, unlike in GR, both shear and 
density gradient tend to zero at early times if 7 > 4/3. 
Thus the high - energy models isotropize into the past for 
realistic equations of state when we include generic linear 
inhomogeneous and anisotropic perturbations. 

IX. CONCLUSION 

In this paper we have given a comprehensive large- 
scale perturbative analysis of flat FRW braneworld mod- 
els with vanishing cosmological constant, extending the 
work presented in [l^ by providing a complete analysis 
of scalar, vector and tensor perturbations for all the im- 
portant stages of the braneworld evolution, namely the 
low energy GR, high energy and dark energy regimes. To 
make this precise we defined dimensionless variables that 
were specially tailored for each regime of interest. In 
this way we were able to clarify further the past asymp- 
totic behaviour of these models and obtain results which 
exactly match the recent work of Coley et. al. 0| 
on the spatially inhomogeneous G2 braneworld models: 
isotropization towards the J-'b model occurs for an equa- 
tion of state parameter 7 > 4/3. 
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